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1 Introduction 



The one dimensional systems, with a time independent Hamiltonian, are treated in a rather different 
way in classical and quantum mechanics. In quantum mechanics the Schrodinger equation for such 
systems determines the spectrum that can include a discrete part for bound states and a continuous 
sector for unbounded (scattering) states. For periodic potentials there are forbidden and allowed 
bands of energy whose borders constitute the spectrum in a wider sense. Only in some cases the 
eigenvalue problem for the stationary Schrodinger equation can be solved "exactly" so that we are able 
to find a "closed expression" for the spectrum (energy eigenvalues) as well as for the corresponding 
eigenf unctions. This situation arises (but it is not the only option) when there are lowering and raising 
operators that close a spectrum generating algebra with the Hamiltonian and factorize the Casimir of 
the algebra up to a constant (BUI 13 Bill]. 

From the point of view of classical mechanics, the Hamiltonians of one dimensional conservative 
systems are integrals of motion, hence all these systems are maximally integrable [6]. This property, 
in classical mechanics, implies that every closed and bounded trajectory must be periodic. For this 
reason, often they are considered almost trivial and usually have not received much attention Q. The 
values of the energy are continuous for both, bound and unbounded motion states. While the motion is 
necessarily periodic with a certain frequency (in general depending on the energy) for bound states, it 
is non-periodic for the unbounded ones. In some special cases, we can find the explicit expression for 
the motion in the phase space by means of the algebraic structure underlying these classical systems. 
The main purpose of this paper is to show how this structure is presented for a class of one dimensional 
classical systems. The Hamiltonians that allow for this treatment are classical analogues of some 
quantum systems. We will see that the algebraic structure of these quantum and classical systems 
are similar, but at the same time there are important differences. Therefore, the general motivation of 
this work is to close the symmetry considerations of both, classical and quantum mechanics in this 
framework. 

The organization of the paper is as follows. In section 2 we introduce our algebraic approach to the 
classical systems in terms of Poisson algebras [8,9, 101- For the factorizable classical systems we find 
two time-dependent integrals of motion similar to the Bohlin invariants lfTTl[T2l[T3l . These invariants 
allow us to obtain algebraically the solutions of motion that give rise to trajectories in the phase space. 
Usually, the time-independent integrals of motion have been used to simplify the differential equation 
of motion in order to get explicit solutions, but the time-dependent integrals of motion have been 
hardly used in studying stationary systems. This work is one step forward in the applications of such 
a type of invariants. This development is systematically applied to some examples in sections 3, 4 
and 5. We will end the paper with some comments and conclusions on the results here obtained. 



2 Factorizations and deformed algebras 

Let us consider the Hamiltonian 

2 

H{x,p) = ^- + V[x) (1) 

2 TO 

where x, p are canonical coordinates, i.e., {p, x} = 1, being {•,•} the notation for the Poisson 
brackets, and V(x) is the potential. Hereafter we will set 2 to = 1 for simplicity. We will investigate 
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a kind of factorization of the Hamiltonian H in terms of two complex-conjugate functions as 

H = A + A- + 7 (H) . (2) 

Here the term j(H), contrary to the usual factorizations in quantum mechanics, may depend on H. 
In the case of bound states with positive energy (H > 0), we will consider A^ linear in p and having 
the form, as suggested from Q, 

A ± = T if{x) p + VH g(x) + <p{x) + 4>(H) (3) 

where the functions f(x), g{x), ip(x) and 4>{H) will be determined in each case. It is clear that 
a global constant factor in the functions A^ will produce an equivalent factorization so that in the 
following we will choose this factor to get simpler expressions. 

We will ask the functions, A^ and H, to define a deformed algebra with the Poisson brackets 
EGO as follows 

{H, A^} = ±i a(H) A ± (4) 

{A+,A-} = -i/3(H). (5) 

The auxiliary functions a(H), (3(H), and <ft(H) will be expressed in terms of the powers of \/H. 
In the case of bound motions with negative energy we must replace the square roots yH by \J — H 
in the above construction. There are systems that can allow for bound and unbounded motions, in 
such cases the complex character of the factors A can change in each energy sector, as well as the 
deformed algebra. This change also happens in the quantum frame for the algebras describing bound 
and scattering states of the same system lfT4l . 

Relation (0]) is the most important since it implies both (O and (f5]). In order to show this, we 
compute 

{H,A + A-} = A + {H,A-} + {H,A + }A- =0 (6) 

where we have made use of (01). From the vanishing of this bracket we conclude that A + A~ will 
depend only on H and therefore, we can express this dependence as in ©. To show ©, we start from 
the Jacobi identity 

{H,{A + , A'}} = -{A + ,{A-,H}} - {A-,{H,A+}} (7) 
substituting on the r.h.s of this identity relation (@]), we find 

{H,{A+,A-}} = -i{A + A-,a(H)} = 

so that the bracket {A + , A~} can be expressed in the form (f5]). 

Now, for a system allowing for this kind of factorization, we can construct two time-dependent 
integrals of motion, 

Q± = A ± eT ia W t . (8) 
It is easy to check that the total time derivative of is equal to zero: 
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since the nonzero partial time derivative is compensated by the poisson bracket of the integral of 
motion with the Hamiltonian by using relation (0]). The (eigen)values of the integrals of motion 
will be denoted 

q± = ce ±id ° (10) 

where c = \q |. Notice that this expression can also change for unbounded motions. Having in mind 
the factorization (0), the modulus c of will depend on the (eigen)value E of the Hamiltonian H: 

c = c(E). 

The two independent integrals Q^(x,p, t) = c(E) e ±td ° allow us to find algebraically the trajec- 
tories (x(t), p(t)) in the phase space. In particular, we can appreciate from ((8]) that the frequencies of 
the motion for bound states are given by a(E). 

If we substitute p 2 = H — V(x) and A given by (O in ((U) we get two relations for the unknown 
functions g(x), f(x), a(H), 4>(H), cp(x) and V(x): 



/(*) = 3^ (^(s) + ¥/(*)) 



-2f'(x)(H - V{x)) + f{x) V'(x) = a(H) ( VH g(x) + <j>(H) + <p(x) 



where the prime denotes the derivative with respect to the corresponding argument of each function. 
If we use A in Q, we have the following equation for "f(H) 

H = f 2 (x)(H - V{x)) + (VHg(x) + <j>(H) + <p(xj) ' + 7 (H) (12) 



which will determine completely the potential. Finally, substituting A^ into (f5]) we will get an equa- 
tion for the remaining function j3{H), 

0(H) = VH[2 f(x) cf{x) - 2 f{x) g(x)] + -j- [fix) g(x) V'{x) + 2 f{x) g(x) V(x)] 

v (13) 
-A fix) <P'iH)(H - V(x)) + 2 f{x) [V'ix) <j>'(H) + <f/(xj\ . 

In the next sections we will study systematically some examples. To get the potentials and the 
functions fix), g(x), (fix), a(H), 7(H) and (3(H), first we have chosen a suitable form of 4>(H), 
g(x) and (f(x), and we have also assumed that H may have a positive or a negative character. 



3 Simple potentials: oscillator, Scarf and Poschl-Teller 

In this section, we will consider the simplest solutions of the problem which are obtained when the last 
term, (f>iH), in the expression © of A^ vanishes. Now, from relation (fl2l we see that there are two 
kind of solutions: those with = 0, ip(x) 7^ 0) and the complementary set (g(x) 7^ 0, <p(x) = 0). 
The first choice leads to the harmonic oscillator, and the second one to the Scarf and Poschl-Teller 
potentials as shown below. 

3.1 The harmonic oscillator ig = 0, if ^ 0) 

This choice gives rise to 

A ± =^if(x)p + ip(x) (14) 
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and using it in (ITTb . we have 

f(x) = ^ r -v'(x) (15) 

- 2 f'(x) (H - V{x)) + f(x) V'(x) = a(H) <p'(x). (16) 

Here, since we assume that H is a variable independent of x, then Eq. ([TBI gives us a(H) = ao and 
Cl5]>-CL6]> take the form 

f(x) = —<p'(x) (17) 

- 2 f'(x) H + 2 f'(x) V(x) + f(x) V'(x) = a ip'(x). (18) 

Taking into account that the coefficient of H in (fT8T ) has to be equal to zero, since it is the only term 
depending on H, from (fTTT ) we have <p(x) = clqx + bo. The remanning part of (fT8T ) gives us 

V(x) = (19) 
Then, substituting y?(a;) in this equation and integrating, we get the harmonic oscillator potential 

V(x) = ^L( a -°f + bo x) + c (20) 
2 a V 2 / 

where c is integration constant. From Eq. (fT2l) . we see that ao = ao/2, and 70 = c — 6q> so t ne 
potential (1201 becomes 

a n 2 



7(x)=^aj + 6oJ +70- (21) 
Using ( fT3l ), we get = «o an< i therefore, in this case we can rewrite all the functions as 

f(x) = l, g(x) = 0, ip(x) = ^fx + b 
a(H) = a , 0(H) = a , -y(H) = 70. 
The explicit expressions of the factor functions are 



(22) 



A ± = T ip+(^x + b Q ) (23) 
and the Hamiltonian is factorized in terms of these functions as 

H = A + A- + 7o . (24) 

Finally, the Poisson brackets read 

{H,A ± } = ±ia A ± , {A + ,A~} = -ia . (25) 

Of course, they constitute the oscillator algebra realized in terms of Poisson brackets. Now, from ([8]) 
we can write the time-dependent integrals of motion in the form 

Q ± = A ± e Tiaot . (26) 
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Having in mind that the Hamiltonian H is an integral of motion whose value is the energy E of the 
system, and the factorization (l24l) . the value of is given by 



q 



: - - yfE^~ Q e ±ie ° (27) 



where the energy is greater than the minimum of the potential, E > 70. Then, taking into account 
(|23T ). (l26l ) and (l27l) . we get the following equations 

Tip+ (yX + 6 ) = y r E^T e ±i{e ° +aot) . (28) 

The phase trajectories (x(t),p(t)) can be found from these equations and have the well known ex- 
pressions 

2 

x{t) = — ( y/E - 70 cos (6»o + a t)- b ) 

ao v y (29) 

= - \/E - 70 sin (6*o + a 1 ) • 

The constant &o can be eliminated by a translation in the x-axis, so hereafter in the next examples it 
will be set equal to zero. The values of fix completely the initial conditions (x(0),p(0)) of the 
motion. 



3.2 The Scarf potential 

Here we will consider (g ^ 0, ip = 0), and the positive character H > of the Hamiltonian for the 
bound states. Following the same procedure as in the previous case, after some calculations we get 
the Scarf potential 

VM= ^m- ■" >0 <30) 

where 70 and a® are constants. We also obtain 

f(x) = cos (^-x), g(x) = sin (^x), cp(x) = 
a(H) = a VH, (3(H) = a VS, 7 (H) = 7o . 
Therefore, we can write the factor functions 



(31) 



= =fi cos (-^-x) p + VW sin {-^-x) (32) 

the factorization, 

H = A+A- + 70 (33) 

and the deformed algebra 

{H,A ± } = ±ia \ r HA ± , {A + , A~} = -ia V~H. (34) 
This algebra can be easily rewritten, defining Aq = VW as follows 

{A ,A ± } = ±i^A ± , {A + ,A~} = -ia A (35) 

which corresponds to the su(l, 1) Poisson algebra. We remark that for the quantum system with Scarf 
potential the spectrum generating algebra is also su(l, 1) Lie algebra 0. 
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Now, the time-dependent integrals of motion satisfying relation (© in terms of are 

Q± = A ± e TiaoVJIt (36) 

with the values g ± = \JE — 70 e ±l 80 . As the energy must be greater than the minimum of the 
potential, E > 70. Substituting (l32l > in (|36"1 ) and taking into account the values of the time-dependent 
integrals of motion, we get 

T i p C o S ^-x + VE sin yi = v 7 ^^ e±*(«o+«o v^t) _ (37) 
From these equations we have the following expressions 



sin —x = y — g — cos (0o + ao v Et) 
p cos ^-x = —^E — jq sin (9q + «o V^E t) 
and the phase trajectories (x(t),p(t)) can be easily written in the form: 



(38) 



2 . / £-7o 
arcsm(Y — ^ — 

y/E(E - 70) sin (0 O + «o v 7 ^ t) 



x(t) = — arcsin(-\/ — — — - cos (6q + ao "v^E t)) 
ao V E 



p(t) 



(39) 



E-(E- 70) cos 2 (0 o + a V#t) 



3.3 The Poschl-Teller potential 

In this section, we assume (g ^ 0, (p = 0), and we choose the negative sign H < of the Hamiltonian 
for the bound states. Therefore, we must change the form of in (O, 

A ± = Tif(x)p + y/^Hg(x) + <p(x) + <f>(H). (40) 

The expressions (fTTb . (fT2l) and ( fT3l ) must also adapt to this situation: replacing the root y/H by \J —H 
in all these equations, and in particular expression (fT3l) takes the form 



(3(H) = V=H [2 /(*) g'(x) - 2 f'(x) g(x)] - -±= [f(x) g(x) V'(x) + 2 /'(*) g(x) V(x)] 
-A f'{x) <t/{H){H - V(x)) + 2 f{x) [V'(x) <//(H) + V '{x)] . 

(41) 

However, later we will also study the unbounded states with positive energy. Following the same 
procedure as in the above cases and having in mind that <fi(H) = 0, we arrive at the Poschl-Teller 
potential 

V{ * ) = -^W*)' ^ >0 (42) 

where 70 and ao are constants. For this case, we get 

f(x) = cosh (^-x), g(x) = smh(^-x), <f{x) = 

(43) 



a{H) = Q V—H, f3{H) = a V—H, j{H) = - 7o 
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Then, we have 

A ± = =Fi cosh ( -^-x) p + V-H sinh (^-x) (44) 

the factorization 

H = A + A- - 70 (45) 

and the Poisson brackets 

{H,A ± } = ±ia V^HA ± , {A + ,A~} = -ia V^H. (46) 



This algebra can also be rewritten, defining Aq = —\/—H as follows 

{A ,A ± } = ±i^A ± , {A+,A-} = ia A (47) 

which corresponds to the su{2) Poisson algebra. Here, we remark that for the quantum Poshl-Teller 
potential the spectrum generating algebra can also be identified with the su{2) Lie algebra |5]. 

Hence, the time-dependent integrals of motion are = A^ e Ttao v~#* s with the values q^ 1 = 
yfE + 70 e ±l e ° . The energy must be negative, but greater than the potential minimum, —70 < E < 0. 
Replacing (1441) in and taking into account the values g , we get 

Ti P cosh yx + V^E sinh % = + 7o e ±i(9o+Q ° v ^* ) . (48) 

Then, we have the following equations 



smh —x = \ — cos (0o + «o y—E t) 

2 V E 

a 



(49) 



p cosh —x = —-\/E — ^q sin (0q + «o y/—Et) 



and from these equations (x(t),p(t)) can be found as 



x(t) = — arcsinh(A / cos (0o + «o V—E t)) 

ot$ \ E 



, yJ-E{E + 10 ) sin (0 O + «o >/ = S t) 
P{t) = / 

-E + (E + 7o ) cos 2 (0 O + a t) 



(50) 



The potential (1421 ) also allows for unbounded motions when E > 0. In this case we can write 
\/—E = i \f\E\ (or \/—H = i y/~H). Then, the factors A^ are no longer complex conjugate, but 
pure imaginary: 

A^ 1 = =pi cosh (—x)p + i v^H sinh (— x). (51) 
The algebra of the Poisson brackets (l46l ) also changes for this energy sector 

{F,i ± } = Ta VHA ± , = a VH. (52) 

But, defining ^4o = — * vH, this algebra has also the form of an su(2) Poisson algebra 

{A ,A ± } = ±i^-A ± , {A + ,A-} = ia A . (53) 
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Dealing with this potential in quantum mechanics, the bound states are described by a finite-dimensional 
unitary representations of su(2) but the scattering states by infinite-dimensional representations of the 
continuous series of su(l, 1) lfl4l . However, if we use the same kind of lowering and raising opera- 
tors, the scattering states belong to infinite-dimensional non-unitary representations of su{2). Here, 
we have seen that in classical mechanics the bound motions are described by the unitary su{2) pois- 
son algebra (in the sense that the factor functions are complex-conjugate) and the unbounded 
motions by non-unitary su(2) algebra (the factor functions are pure imaginary) by means of the 
Poisson brackets. 

The integrals of motion are imaginary and their eigenvalues can be written in the form q^ 1 = 
=pi y/E + 70 e ° . After similar calculations, these integrals of motion have the explicit expressions 

TP cosh yx + y/E sinh ^-x = + 7o e^ e ° +a ° (54) 



Therefore, we get 



sinh ^-x = J — sinh (9 + ao \fE t) 



2 

Finally, the unbounded non-periodic solutions are 



(55) 



p cosh ^-x = \JE + 70 cosh (6q + ao VEt) 



2 . h( E+jo 

— arcsmhl \ / — — 

*o V E 

y/E(E + to) cosh (6o +a VEt) 



x(t) = — arcsmh(-\/ — — — - sinh (^o + ao ^E t)) 
a V E 



Pit) 



(56) 



E + (E + 70 ) sinh^(#o + a VEt) 



Remark the different dependence of the frequencies on energy for the bound motion states of the 
above three systems. In the oscillator case, it is a constant independent of E, for the Scarf potential 
the frequency is an increasing function of E, while in the Poschl-Teller potential the frequency is a 
decreasing function of E. This corresponds to the values of the bracket (@]) of the Poisson algebra for 
each of these systems. The potentials and some trajectories are plotted in Figures [HE] and [3] We also 
notice that the values of 70 in the potentials (l30l) and (|42l> have been taken with suitable signs giving 
rise to bound states. 



V(x) 




Figure 1: The oscillator (continuous line), Scarf (dashed) and Poschl-Teller (dotted) potentials. 
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x(t) 



Figure 2: Some phase trajectories for the same three energies of motions in the oscillator (a) and Scarf 
potentials (b). The exterior trajectories of the latter have higher frequencies than the inner ones. 
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Figure 3: Now, it is shown the phase trajectories for the same three energies of periodic motion in the 
Poschl-Teller potential. The exterior trajectories have slower frequencies than the inner ones. Also it 
is shown some unbounded trajectories. 



4 Singular potentials 



In this section we consider <p(H) / 0, to obtain other classical systems. Here, we also propose some 
specific form of <fi{H) and consider j(H) different from constant so that we will find the potentials 
of the previous section including a singular term. 



4.1 The singular oscillator 

If we choose <f>(H) = (3q H, and (g(x) = 0, tp(x) / 0), after similar calculations to the preceding 
cases, we get the singular oscillator potential 

v ^ = ((3 a 2 x-4b ) 2 + 4ag( 7 o + dg) _ do (5?) 



16%a 2 (Poa'x-Uo) 2 ft 



o 
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where (3 , a > an d do are constants. Henceforth, we will choose the global constant f3 = 1, the 
potential origin do = 1/2, and set the space origin 60 = 0, so that we will handle the simplified 
expression 

V(x) = ^f + ^L±l-l (58) 
lb Og x z 2 



For this choice, the unknown functions take the form 



f(x) = -\alx, if{x) = -^ + \ 



13(H) = 2a H - a , 7(H) = —H 2 + 7o 
with a(H) = ao- Hence, the explicit expressions for the factors are 



(59) 



.4. aox cxkx 2 1 

A ± = ±i^-p--^ + H+- (60) 
and the relation between the Hamiltonian and these functions is given by 

H = A + A- -H 2 + 7o . (61) 
Finally, we have the su(l, 1) Poisson algebra 

{H,A ± } = ±ia A ± , {A + ,A~} = -ia (2H -I). (62) 



In this case, the time-dependent integrals of motion has the form d26l ). with c(E) = \J E 2 + E 
Then, the phase trajectories (x(t),p(t)) for E > ~ 1 +v /1 + 4 t o are 



7o- 



x(t) = — y/2E + 1-2 c{E) cos (6» + a t) 
ao 

(63) 

_ c(E) sin (gg + q t) 

P ~ + 1-2 c(£) cos(fl + a 1) ' 

These expressions show that the frequency of the motion for the singular oscillator is a constant ao 
independent of the energy E. 



4.2 The generalized Scarf potential (H > 0) 

For the specific form 



<P(H) = ^=+(3 Vh (64) 

V ti 



and tp(x) = 0, we have the generalized Scarf potential. Here we will assume @o / 0, and by means 
of a global constant of the factor functions we can set it equal to one: /3q = 1. In this way we get the 
potential 

W "4Uin 2 (^x) cos 2 (^x) (Mj 
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where 70, ao> an d So w& constants. We have also 

f( x ) = -2 cos(^x) sin(^x), g(x) = -2 sin 2 (^-x), ip(x) = 

A 2 A 2 (66) 

P(H) = a y/H--^%, j(H) = 70 - % a(H) = a VH. 

Then, the two complex conjugate functions are 

A ± = T2i cos (— x) sin (— x)p + 2 V# sin 2 (^-x) - -% - VH (67) 
M 4 4 



with the factorization relation 

# = + 70 - | (68) 

and brackets 

{i7,A ± } = ±ia Vl : ^ ± , {A + ,A~} = -ia (^y/H - j^=^j . (69) 

In this case, as well as in the remaining examples, the Poisson algebra is not a Lie algebra, but a 
deformed algebra due to the nonlinear function in the last bracket of d69"l ). We can write the time- 
dependent integrals of motion in the form = e ^ ta oVHt^ w j t j 1 ^ va i ue g± = c(E) e ±l9 ° 
where 

c(E) = y_ 7o + J E + |L. (70) 

In order to have a positive potential with bound states, if we assume 70 > 0, then the coefficient So 
must satisfy 

-7o<25 <7o- (71) 

The energy of the bound states are given by the constraint c(E) > 0, that is, when E is bigger than 
the minimum of V(x), 

E> . 

Then, the trajectories (x(t),p(t)) are periodic for any allowed value of E: 

2 ( S + c(E)VE cos (9 + a VEt) , 

x{t) = — arccos 



ao E 

(72) 



c(E) E sin (6»o + a V E t 

p(t) = - 



E 2 - (So + c(E) VE cos (0 O + «o y/E t)) 2 



4.3 The generalized Poschl-Teller potential (H < 0) 

In this case, we consider H < and accordingly, we change the sign of H in the square roots of the 
function 4>(H) = + Po V~H, keeping (p = 0. We also make use of expression (l40l of A^, and 
equation (|4TT) of Section 3.3. In this way we get the generalized Poschl-Teller potential 



V< ^ - ± I 2 ^° + 70 + 2 ^° ~ 70 I (IV* 
{X) < 1 sinh 2 {^-x) cosh 2 (^x) 1 
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(74) 



and the functions 

f(x) = 2 cosh(-^x) sinh(-^x), g(x) = 2 sinh 2 (-^-x), f(x) = 

(3(H) = -a V^H- "° ° , 7(fl-) = 7o + 2 J ff + ^, a(H) = a V^H 

where ao, So, 7o are constants, and we have set (3q = 1. Now, the factors A 1 * 1 are given by 

A± = ^2 i cosh (^x) sinh (^x)p + 2 V^tf sinh 2 (^x) + -=L= + (75) 
4 4 4 y — 

with the factorization relation 

-H = A+A- + 7o + § (76) 
and the deformed Poisson algebra takes the form 

{H,A ± } = ±ia V^HA ± , {A + ,A~} = ia (—^j== + \ . (77) 

The corresponding integrals of motion in terms of A^ 1 are = e Tiao ^~ Ht j with = 
c(E)e ±i9 °, where 

<E) = ^-70 - E - | . (78) 

If the potential has a positive singularity at x = 0, in order to allow for bound states, the coefficients 
in the potential, besides (TTTb . they must satisfy 1 2 <5o + To I < |2#o — 7o|- Then, the negative energy 
values of such states are given by the positive character in the square root (|78T > of c(E), 



^ + V^-M <E<0 . (79 ) 



2 

Then, for these energies, the periodic trajectories (x(t),p(t)) have the following form: 

ao E ' 

c(E) E sin (6»o + a \f zr E t) 

p(t) = - 



(80) 



(S - c{E) V^Ecos (0 O + a V^Et)) 2 - E 2 



In order to get solutions for the unbounded motion, we choose E > 0. Then, the integrals of 

motion are imaginary and the eigenvalues can be written as = i c(E) e T e ° . In this case 

the algebra (1771 ) changes its character. After similar calculations, we have the unbounded phase 
trajectories (x(t),p(t)): 

x n\ = 2 arccoshf fo + ^ cosh ( g ° 
ao E 

c(E)E smh(eo + a VEt) (81) 



Pit) 



(So + c(E) y/E cosh (6 Q + a VE t)) 2 - E 2 
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where 



(82) 



The figures of these singular potentials as well as some trajectories in the phase space are depicted 
in Figures 0] and [51 Notice that the frequency of the periodic motions in the above potentials is again 
governed by a(E), and it has the same expression as in the non-singular potentials of section 3. 



V(x) 



V(x) 



V(x) 




\ ; 

\ <' 10 


1 1 




V ' 


y 




\* ' 






\ 5 




b) 



Figure 4: Plot of the oscillator (a), Scarf (b) and Poschl-Teller (c) potentials (continuous lines) to- 
gether with their singular counterparts: the singular oscillator, generalized Scarf and generalized 
Poschl-Teller potentials. 




Figure 5: Plot of three closed and three open trajectories in the phase space for the generalized Poschl- 
Teller potential. 



5 One dimensional Morse potentials 

In order to get the one dimensional Morse potential and other related potentials, we will consider in 
this section the special form 4>{H) = -4fi= with Lp = 0. This choice corresponds to the case excluded 
in the previous section, when (3q — > 0. 
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5.1 Hyperbolic case (H < 0) 

For <p(H) = J^jj , taking into account (l40l > and (|4TT) . we get the potential 

tti , 2J (Ce> + De->)-70 

V(x) = S5 (83) 

{Ce**-De-^ x y 

where ao, do, 70 and C, D are constants. The functions f(x), g(x), (3(H) and j(H) have the 
following form 



f(x) = Ce^ x - De~^ x , g(x) = C e^ x + D e~^ x , <p(x) = Q 



(84) 



with = «o V—H. Then, the poisson brackets are 



ao 5n 



{H,A ± } = ±ia y/=HA ± , {A + ,A-} = i( — ^= + 4, C D a V^H) . (85) 

The functions .4. can be written as 

A ± = Ti{Ce^ x -De~^ x )p + V^H{Ce^ x + De-^ X ) + -^L= (86) 

\/-H 

and the time-dependent integrals of motion 

q± = A ± e Tiao ^ 1It , q ± =c{E)e ±ie °. (87) 

We shall consider different cases according to the values of C and D, leading to regular potentials 
with bound states. 

5.1.1 C = -D = 1/2. 

The option C = — D is equivalent to any other with opposite signs of C and D; it is enough to 
consider an x-translation. By means of a constant factor we can choose C = — D = 1/2 then, the 
corresponding potential is known as Scarf II (hyperbolic) @ 

V(x) = 2 5 tanh(^x)sech(^x) - 7o sech 2 (^x). (88) 

We have the following relations obtained for these particular values of C, D: 

H = A+A- - 70 + f (89) 

{H,A ± } = ±ia V^HA ± , {A + , A~] = i a ( tt f== - (90) 

H \/—H 

For this case the functions A^ read 

A ± = Ti cosh(^x) p + V^HsmhA) + -=L= . (91) 
2 2 \J — H 
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Then, (x(t),p(t)) have the following form: 



and 



p{t) = _ c{E)E sin (e + a V=Et) ^ 

'e 2 + [5 - c(E) V^Ecos (0 O + a V^Etrf 



s 2 

where c(E) = a/ 70 + E — -Sr. From c(E), we see that the range of the energy for bound states is 



given by (70 — yTq + 4<5q)/2 < i? < 0. When i? > 0, we have unbounded motion states with pure 
imaginary functions Ar and a different algebra from d90T >. similar to that of Poschl-Teller potential. 

For the choice C = D, we get a trigonometric singular potential without bound states. 
5.1.2 D = 0. 

In this case, by taking D = 0, the corresponding potential and the functions /(x), <?(x), (3(H), and 
7(i?) can be obtained from (f83T > and d84l ). We can also set C = 1, since C is a global non-vanishing 
constant factor. The result is the Morse potential 

V{x) = 2 5 e-^ x --f e~ aox . (94) 



Then, the complex conjugate functions are 



A ± = ^ie^ x p + V^He^ x + -^=. (95) 

y/-B 

Now, we can build the usual time-dependent integrals of motion with eigenvalues = c(E) e^ 1 1 e ° 

n? — 

where c(E) = y — + 70. The periodic trajectories (x[t),p(t)) have the following form: 

2 r 5 -c(S)^/^Ecos(0 o + ao >/ = St) n m ^ 

x(t) = — log (96) 

ao E 



and 



(f) = c(E)E sin (9 + a V=Et) 

5 - c(E) V^E cos (0 O + a V^E t) ' 



In order to have periodic bounded motions we must take the coefficients in the potential with the 
signs 5o < and 70 < 0. In this way the range of the energies is 5 2 /70 < E < 0. 

The case C = gives similar results, where the exponentials have opposite signs than in the case 
just considered above. 

5.2 The trigonometric case (H > 0) 

When we take 4>(H) = 4|=, we obtain the potential 

V ix) = - ^ C 't" + De ~V +7 ° (98) 
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V(x) 





Figure 6: Plot of the hyperbolic Scarf II potential with the parameters 70 
Morse potential with 70 = —4, 5q = —4, on the right (b). 



4, 5 = 2 (a), and the 



where ao, do, 70 and C, D are constants. The involved functions have the following expressions 

g (x) =Ce i ^ x + De~ i ^ x , <p(x) = 



f(x) =i(Ce l 2 x -De 
0(H)- 



-^L + 4CDa VH, 7(ff)=7 + (l-4CL>)#-f 



(99) 



with = aovH. Then, the associated Poisson brackets are 

{ff,i ± } = iia \/Ii ± ) {A + ,A~}=i( 
For this case take the form 



«o^o 



ACDao^H). (100) 



1 ■ a o ■ a o 1 — ■ a o ■ a o on 

A ± = ±(Ce l ^ x -De- l ^ x )p + VW(Ce l ^ x + De- l ^ x )+ " 



H 



(101) 



and the time-dependent integrals of motion are given by = e T 1 a " v^*, with = c(E) ' 
When we take D = —C = i/2, we obtain the potential called Scarf I (trigonometric) 



V(x) = 25 tan (^-x) sec(^x) + 70 sec(^x) 



(102) 



together with 



f(x) = cos(^x), g(x) = sin(fix), <p(x) = 



Then, the Poisson brackets are 



y + aoVH, 7(-ff)=7o-f, a{H) = a VH. 



(103) 



Now the functions yl^ take the form 



A = =F« cos(— x) p + V ri sm(— x) H — = 
2 ^ \f H 



(104) 



(105) 
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and 



c(E) 



E + 



E 



70 • 



(106) 



From this expression we rind that bound motion states exist when 70 > 25q. Under this condition, the 
allowed energies are E > (70 + a/To" 
difficulty as in the Scarf potential. 



4<5q)/2. The corresponding motions can be found without 




Figure 7: Plot of the trigonometric Scarf I potential for the values: <5 = 1, 7o = 4. 



6 Conclusions 

In this work we have studied a whole class of one-dimensional classical systems characterized by an 
underlying Poisson algebra which in general is a deformed Lie algebra. Here, the Poisson algebra is 
not made up of (time-independent) integrals of motion, as it is the usual case ['8,9, 10 ], but it includes 
functions, A^, directly related with time-dependent integrals of motion. 

We have obtained some systems like the Scarf, Poschl-Teller, Morse etc. which are clearly the 
classical analogues of the one-dimensional quantum systems that can be solved by means of the 
factorization method. However, the factorizations here employed are not exactly the factorizations 
of the corresponding quantum cases [2]. Instead, they are more related with the so called spectrum 
generating algebras in quantum mechanics lfl"4l l3ll5Tl. as can be seen replacing Poisson brackets by 
commutators: {•, •} — > ■]. Therefore, we have shown that these algebras can also be useful in 
classical mechanics to compute time-dependent integrals of motion of Bohlin type that give us the 
solutions of the motion in the phase space lfTTl[T2l[T3l . 

The algebraic structures here obtained for some classical systems correspond to well known ones 
for the analog quantum mechanical systems. Hence, this correspondence is very important to describe 
the coherent states of such quantum systems: for instance, the expected values of x and p are adequate 
for the harmonic oscillator, but in the case of the Scarf potential we should consider the expected 
values of sin x and p cos x, which are the components of the functions , as we see from d29l and 
OH), respectively B1ID1. 
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